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ABSTRACT 


Th«  performance  characteristics  of  a  residential  eolar 
assisted  water-to-air  heat  pump  heating  system  are  analysed 
for  State  College,  Pennsylvania .  A  realistic  residence  and 
solar  assisted  vater-to~alr  heat  pump  system  are  modeled  for 
this  northern  climate  using  the  transient  simulation 
computer  code  TRNSYS  developed  by  the  University  of 
Wisconsin*  The  system  Is  studied  over  a  one  month  winter 
period,  December,  using  actual  hourly  weather  data  for  State 
College*  The  eystem  is  analyzed  for  both  the  cloudiest  and 
clearest  December  weather  recorded  in  the  last  30  years. 
The  collector  area  and  storage  tank  capacity  are  varied  and 
the  effects  on  system  performance  are  studied.  Other 
parameters  verled  and  studied  are  the  collector/  storage 
tank  heat  exchanger  effectiveness  and  the  room  temperature* 
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CHAPTER  l 


1.1  Introduction 

With  the  energy  erlsie  vorsening  end  eonvencionel 
heeclng  fuels  becoming  scarce  and  expensive  a  need  exists 
for  the  development  of  alternative  heating  fuels  and 
systems.  Solar  energy  which  Is  "free"  and  readily  available 
la  a  promising  source  of  building  heat.  Unfortunately  for 
northern  climates  direct  solar  heating  of  buildings  can 
provide  only  a  small  portion  of  the  total  heating 

requirements  without  extremely  large  collector  areas  and 
heat  storage  capacity ( 1  ]  •*  The  solar  energy  can  however  be 
used  to  augment  conventional  heating  systems  or  used  to 
Improve  the  performance  of  the  system,  as  with  a  heat  pump. 
The  heating  capacity  of  the  heat  pump  decreases  as  the 
evaporator  temperature  <l.e.  the  outdoor  ambient  air 

temperature)  decreases.  At  the  same  time  the  heating  load 
of  the  residence  Increases.  The  point  at  which  the  heat 
pump  capacity  equals  the  heating  load  Is  called  the  balance 
point.  The  heat  pump  Is  an  economical  heat  source  as  long 
as  it  can  operate  at  a  temperature  above  Its  bslance  point. 
In  cold  northern  climates  the  evaporator  inlet  temperature 
of  the  standard  air-to-air  heat  pump  installation  Is  often 
below  the  balance  point.  For  water-to-air  heat  pumps  an 


*  Numbers  in  brackets  refer  to  references  listed  at  the  end 
of  this  paper. 
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existing  warn  water  source  (teaparature  above  the  balance 
point)  for  the  heat  puap  evaporator  la  difficult  to  find  In 
northern  cllaates  particularly  for  residential 
Installations*  Solar  energy  can  be  uaad  easily  to  provide 
this  wara  taaparatura  source  for  the  hast  puap  evaporator 
sines  only  a  low  grads  teaparature  (50-100  degress  P)  is 
required . 

Over  the  past  several  years  aany  papers  and  studies 
have  been  done  on  the  evaluation  of  solar  energy  heat  puap 
systeas .  The  earliest  research  on  solar  assisted  heat  pumps 
was  begun  by  the  American  Gas  and  Electric  Service 
Corporation  In  1949  In  a  5  year  study  which  used  a  high 
temperature  storage  system  with  the  heat  puap  coupled 
directly  to  the  collector  1 2) .  In  1952  R.  C.  Jordan  and  S. 
L.  Threlkeld  published  a  paper  in  which  the  potential  of 
this  system  on  the  basis  of  some  experiments  was 
described ( 3] .  Since  that  tine  the  solar  assisted  heat  pump 
system  has  been  improved  upon  and  extensively  studied.  The 
effects  of  geographical  location  and  variations  in  the 
system  design  or  component  design  have  bean  addressed  In  the 
present  literature [4-6] . 

The  study  described  In  this  paper  will  analyse  the 
performance  of  a  solar  assisted  vater-to-alr  heat  puap 
residential  heating  system  in  State  College,  Pennsylvania 
(latitude  40.5  degrees  N).  A  realistic  realdencc  and  solar 

i 
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assisted  vattr-to-alr  heat  pump  system  ere  modeled  to  study 
the  performenee  characteristics  o£  such  a  system  in  e 
northern  climate.  The  transient  simulation  computer  code 
TRNSYS  developed  by  the  University  of  Wisconsin  and  as 
modified  at  the  Pennsylvania  State  University  is  used.  The 
system  is  studied  over  a  one  month  winter  period,  December, 
using  actual  hourly  weather  data  for  State  College.  The 
system  la  analyzed  for  both  the  cloudiest  and  clearest 
December  weather  recorded  In  the  last  30  years.  The 
collector  area  and  storage  tank  capacity  are  varied  and  the 
effects  on  system  performance  are  studied.  Other  parameters 
varied  and  atudlad  are  the  collector/storage  tank  heat 
exchanger  effectiveness  and  tha  room  temperature. 


CHAPTER  2 


2.1  Description  of  Residential  Structure 

For  investigation  end  Modeling  purposes  the  solsr 
assisted  heat  pump  system  was  assumed  to  be  installed  in  a 
single  family  residence.  It  is  a  two  story  structure  with  a 
full  basement  and  attached  two  car  garage.  The  living  area 
is  approximately  2200  square  feet  with  8  foot  ceilings 
throughout.  The  solar  collectors  are  located  on  the  south 
slope  of  the  roof  and  are  tilted  at  50.5  degrees  (latitude  * 
10  degrees),  Che  recomended  [2]  optimum  angle  for  heating 
applications.  A  floor  plan  of  the  residence  is  shown  in 
Figure  1.  The  residence  is  located  in  State  College, 
Pennsylvania, ' a  small  university  oriented  city  situated  in 
the  center  of  the  state. 

An  energy  conserving  design  has  been  utilized  to 
minimize  Che  required  heating  load.  The  house  walls  are 
built  of  wood  franc  construction  with  four  inch  face  brick* 
All  windows  are  double  pane  glass  and  a  stora  door  is 
Installed  on  the  front  door*  Three  and  one  half  Inches  of 
1-1 1  insulation  is  installed  in  the  walls  and  six  Inches  of 
B-22  insulation  is  placed  la  the  second  floor  celling.  A 
aore  complete  list  of  the  construction  details  can  be  found 
In  Appendix  A* 
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In  order  to  also  and  aalact  tha  appropriate  heat  puap 
It  vai  necessary  to  calculate  a  daalgn  auaaar  cooling  load 
aa  vail  aa  tha  winter  heating  load.  Tha  heating  load  vaa 
calculated  uaing  tha  general  proeaduraa  aat  forth  by 
ASHRAE(31  uaing  ASHRAE  tabulated  data.  Tha  resulting  daaign 
heat  loaa  vaa  approximately  32222  Btu/hr.  Fro*  this  tha 
building  UA  (tharnal  tranaalt tanca  X  exposed  area)  vaa  found 
to  be  311  Btu/hr-F  which  indicataa  a  vary  wall  Inaulated 
houaa.  Tha  cooling  load  vaa  dateralnad  uaing  the  ASHRAE 
slapllfled  procedures  for  residential  calculations [4]  and 
ASHRAE  tabulated  dat*.  Tha  cooling  load  requirement  was 
found  to  be  approximately  1.4  tons.  Suaaary  calculation 
sheets  for  both  tha  heating  and  cooling  loads  can  be  found 
in  Appendix  A. 

2.2  System  Description 

A  schematic  diagram  of  the  residential  solar  assisted 
hast  pump  heating  system  is  shown  in  Figure  2.  The  major 
components  of  the  system  are  a  flat  plate  solar  collector, 
storage  tank  and  a  vater-to-alr  heat  pump.  Other  equipment 
Includes  pumps,  heat  exchanger,  a  relief  valve  and  auxiliary 
electric  heating  which  is  an  Integral  part  of  the  heat  pump. 
The  major  components  are  discussed  in  more  detail  la  the 
following  paragraphs. 


FIGURE  2  -  Schematic  of  solar  assisted 
hast  pump  hatting  system 


A*  Flat  Plata  Collector  -  The  flat  plata  aolar  collector 
chosen  for  tha  simulation  van  a  Pittsburg  Plata  Glasa(PPG) 
Model  CD414.  The  PPG  collactor  features  weather  sealed 
construction  with  two  tempered  glass  cover  plates  and  an  all 
copper  absorber  plata  with  a  non-select lve  coating*  Tha 
collactor  la  especially  designed  for  applications  where  tha 
collactor  area  la  Halted  or  where  large  temperature 
differentials  exist  between  the  aablent  air  end  the  absorber 
plate*  Each  collector  panel  Is  approximately  20  square 
feet*  The  collector  array  uses  a  40%  ethylene  glycol-water 
solution  as  the  circulated  heat  transfer  aedlua*  This 
provides  protection  froa  freeilng  to  -12  degreee  F* 
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ADDITIONAL  PARAMETERS 

Wind  Velocity:  0-5  mph 

Fluid:  401  ethylene  glycol  solution 

Cp  -  .792  Btu/lb-F 

Flowrate  -  .025  gpa/sqft 

Tube  Diameter  •  .5  In 

Tube  Spacing  ■  3  in 


•X  #2  *3  »5  *6  *7  *9  X*0 

(  Tf 1  -  T»  )  /  HT 


FIGURE  3  -  Efficiency  of  the  PPG  flat  plate 

aolar  collector 


The  perforaance  of  this  flat  plats  collector  can  be 
described  by  a  nuabcr  of  paraaeters.  Perhaps  the  aopt  basic 
perforaanee  paraaeter  Is  ehe  collector  efflcleney(  )  which 
for  the  PPG  collector  la 


-.69-.7(<Tf 1-Ta)/HTJ  (1) 

where  Tfi"lnlet  fluid  temperature (degrees  F) 
Ta  "ambient  air  temperature (degrees  7) 
NT  "incident  radiation  per  unit  area  on 
the  collector  aur f ace(Rtu/hr-sqf t) 


The  efficiency  equation  la  shown  graphically  In  Figure 
3  along  with  a  number  of  other  parameters  describing  the 
construction  of  the  PPG  flat  plate  collector*  These 
parameters  Include  the  transmittance  (?)  of  the  cover,  the 
abaorptance(et)  of  the  absorber  plate  and  the  emmlsslvlty (£) 
of  the  absorber  plate* 

Three  other  useful  parameters  describing  the  PPG 
collector  are  the  Collector  Efficiency  Factor(FR),  the 
Collector  Overall  Energy  Loss  Coef f lelent (UL)  and  the 
Collector  Geometry  Efficiency  Factor(F').  These  perameters 
are  related  to  the  Rate  of  Useful  Energy  Collected  per  Unit 
Area  of  Collector (qu)  according  to  the  following  expression 
developed  by  Hottel  and  Woerts(5]  and  extended  by 
Uhllller [6]  and  Bllss(7]. 

qu-FR(HT<rK)-UL<Tfi-Ta)l  (2) 

Since  aqu/HT  equatlon(2)  can  be  rewritten  as 

-Ft(?o<)-tFR*UL(Tf  i-Ts)/HTJ  (3) 

Comparing  this  to  ths  collector  efficiency  equation  for 
the  PPG  collector.  It  follows  that 


FR(?ot)-*69  and  FR*UL-.7 


(4) 


Substituting  known  valuta  for  tho  transmit  tones  and 


absorptance 


FR-.69/(.72) <.95)-1.0 

thua 

UL-.7  Btu/hr-oqf t-F 


(5) 

(6) 


Finally,  using  UL,  the  Collector  Geonetry  Factor(F*) 
can  be  obtained  fron  graphs  developed  by  Duffle  and 
Becknan  [8] •  The  value  of  F '  for  the  PPG  collector  is  .99* 

B.  Storage  Tank  -  An  insulated  steel  underground  tank 
containing  water  is  used  for  theraal  storage.  The  overall 
theraal  transalt tance (0 )  of  the  tank  la  .05  Btu/hr-aqf t-F. 
The  tank  la  cylindrical  in  shape  with  a  constant  length  to 
dlaaeter  ratio  of  three(as  the  voluae  is  varied  in  the 
siaulatlon).  A  tube-ln-she 1 1  heat  exchanger  with  an  assuaed 
effectiveness  of  70S  is  used  to  transfer  heat  froa  the 
collector  loop  to  the  storage  tank. 

C.  Water-to-Alr  Heat  Puap  -  An  Aaerlcan  Air  Fllter(AAF) 
water-to-alr  heat  puap  Model  HW19  (1*5  tons)  was  selected 
for  the  residence.  As  is  standard  practice  the  heat  puap 
was  selected  to  aeet  the  cooling  load  requlreaenta  of  tha 
building.  It  is  recognised  that  in  harsh  weather  the  heat 
puap  will  be  unable  to  provide  100X  of  the  heating  load 
therefore  the  AAF  heat  puap  is  equlppad  with  a  parallel 


auxiliary  electric  raalatanca  haatar  Co  provlda  tha  deficit. 
Tha  AA F  Modal  HW19  vater-to-air  haat  pump  la  raced  Co 
provide  600  CFM  of  heaced  air  co  cha  house.  Tha  water  flow 
on  cha  water  aide  of  Cha  haaC  pump  la  5  gpa.  Table  1  la  Cha 
haaclag  performance  charac carls t ica  of  cha  Modal  HV19  for 
aalacalalag  Che  room  cemperacure  at  70  degrees  P* 


TABLE 

1  -  Heating 

performance 

characterletlcs 

of  AAF  Model  HW19 

has  t 

pump 

EVAP 

HEAT 

HEAT 

INLET  TEMP 

ABSORBED 

REJECTED 

POWER  INPUT 

COP 

(degrees  F) 

(MBtu/hr ) 

(MBtu/hr) 

(kw) 

(MBtu/hr ) 

60 

11*3 

18.0 

1.80 

6.16 

2.93 

30 

13.3 

20.7 

2.06 

7.03 

2.96 

60 

13.2 

23.2 

2.30 

7.85 

2.96 

70 

17.  1 

23.9 

2.33 

8.63 

3.00 

80 

19.0 

28.6 

2.76 

9.62 

3.01 

90 

20.9 

31.1 

2.99 

10.20 

3.05 

100 

22.8 

33.7 

3.22 

10.99 

3.07 

110 

26.7 

36.3 

3.63 

11.77 

3.08 

2*3  Operation  of  System 


Tha  solar  asalatad  hast  pump  system  la  designed  Co 
supply  cha  complete  heating  raqulramants  of  cha  raaldence 
with  minimal  energy  consumption*  Tha  required  haac  is 
provldad  by  one  of  three  modes*  Thesa  modes  are  direct 
haaclng  from  storage,  haac  pump  heating  and  auxiliary 
electric  resistance  heating. 
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The  collector  loop  usee  a  402  ethylene  glycol-water 
mixture  to  prevent  freexlng  during  low  ambient  temperature 
conditions*  The  mixture  Is  pumped  through  the  collector/ 
heat  exchanger  loop  at  a  rate  of  *025  gpm/sqft  of  collector* 
The  pump  Is  controlled  by  a  control  mechanism  which  turns 
the  pump  off  when  the  temperature  of  the  fluid  leaving  the 
collector  la  less  than  the  storage  tank  temperature.  The 
pump  resumes  operation  when  the  temperature  of  the  fluid  in 
the  collector  reaches  the  storage  tank  temperature  plus  10 
degrees  F.  This  control  sequence  prevents  the  loss  of 
thermal  energy  from  the  tank  by  reradlatlon  through  the 
collector  at  night  or  cloudy  days  when  solar  insolation 
decreases  to  a  minimum* 

When  the  house  thermostat  calls  for  heat  and  the 
storage  tank  temperature  is  greater  than  or  equal  to  110 
degrees  F  water  from  the  tank  is  pumped  through  a  heating 
coll  in  the  house  air  duct  while  the  house  air  fan 
circulates  warm  air  through  the  house*  This  is  direct 
heating  from  the  storage  tank.  If  heating  is  required  and 
the  storage  tank  temperature  is  below  110  degrees  F  the 
direct  heating  mode  is  disabled  and  the  control  valves 
circulate  the  water  through  the  heat  pump  heat  exchanger* 
The  heat  pump  now  provides  heat  to  the  house*  On  occasion 
the  capacity  of  the  heat  pump  will  be  unable  to  meet  the 
house  heating  requirement  In  which  case  the  built-in  house 
air  stream  electric  resistance  heater  of  the  heat  pump  will 


provide  the  additional  required  heat*  If  the  storage  tank 
temperature  fella  below  40  degreea  F  end  heating  le  atlll 
required  then  the  heat  pump  is  shut  down  and  all  the 
required  heating  Is  provided  by  the  electric  resistance 
heater. 

If  the  aablent  air  temperature  exceeds  65  degrees  F  the 
heating  requirements  decrease  to  zero  and  all  heating  modes 
become  Inoperative. 

2.4  Computer  Model 

The  aolar  assisted  heat  pump  system  shown  In  Figure  2 
was  modeled  using  the  Transient  Simulation  Program  (TRNSYS) 
developed  at  the  University  of  Wisconsin.  The  program 
allows  one  to  make  a  transient  (time  dependent)  simulation 
of  the  overall  system  by  interconnection  of  Individual 
component  models.  The  TRNSYS  modeling  components  used  sre 
shown  In  Figure  4  which  Is  e  flow  diagram  of  the  computer 
model . 


The  reader  component  Is  used  to  input  hourly  weather 
data.  The  simulation  Is  run  for  one  month  of  the  heating 
season.  Only  ambient  temperature  and  solar  Insolation 
(measured  on  a  horizontal  surface)  data  are  needed  for  the 
simulation.  Actual  weather  data  obtained  from  the 


Pennaylvanla  State  University  weather  station  for  the  months 


of  December  1960  and  December  1972  an  used.  For  Che  past 
thirty  years  of  available  weather  data,  these  years 
represent  the  highest  and  the  lowest  values  for  December 
solar  insolation,  respectively,  for  State  College, Pa. 

The  radiation  processor  component  converts  the  solar 
insolation  or  radiation  on  a  horixontal  surface  (H) ,  to 
solar  radiation  on  a  tilted  surface  (HT).  The  collector  la 
tilted  at  an  angle  of  50.5  degrees.  By  using  Mode  3  of  the 
radiation  processor  component  the  total  radiation  on  a 
tilted  surface  la  a  combination  of  the  beam  radiation, 
diffuse  radiation  and  reflected  radiation.  The  radiation 
processor  calculates  these  values  from  the  horizontal 
surface  radiation,  the  collector  orientation  and  the  ground 
reflectance  which  la  .2  for  bare  ground. 

The  liquid  collector-storage  module  la  a  single 
component  used  to  model  the  solar  collector,  ethylene 
glycol/water  heat  exchanger,  storage  tank  and  associated 
pumps  and  relief  valve.  A  schematic  of  the  module  Is  shown 
In  Figure  3.  The  combining  of  these  components  Into  one 
subsystem  reduces  computation  time  and  expense.  The 
collector  fluid,  40Z  ethylene  glycol  solution,  is  pumped  at 
the  collector  manufacture rs  recommended  rate  of  .023 
gpm/sqft  of  collector.  Across  the  heat  exchanger  the 
capacity  rate  ratio  is  held  at  1.0  [i.e.  (mCp)eg  -  (mCp)v  ]. 
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The  Initial  storage  tank  temperature  at  the  beginning  of 
December  la  assumed  to  be  70  degrees  F.  It  turns  out  that 
any  Initial  temperature  in  this  range  gives  nearly  identical 
results  over  the  one  month  period*  The  relief  valve  Is  set 
to  maintain  the  tank  temperature  belov  212  degrees  F* 


FIGURE  5  -  TRNSYS  liquid  co 1 lec t or / s tor age  module 


The  water-to-alr  heat  pump  component  uses  actual  heat 
pump  performance  characteristics*  The  model  Is  quasi-steady 
state  In  nature.  That  is,  at  any  Instant  in  the  "transient" 
simulation,  the  heat  pump  performs  as.  if  it  were  at  steady 
state,  having  for  the  moment  constant  source  and  sink 
temperstures  and  hence  a  unique  heat  addition,  heat 
rejection  and  work  Input*  The  heat  pump  Is  modeled  to 
operate  In  the  heating  mode  only.  A  parallel  auxiliary 
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The  quantity  Integrator  almply  keepa  a  eontlnuoua  total 
of  varloua  ayaten  outputa  throughout  the  simulation*  It 
vorka  alolllar  to  a  KWH  meter  which  contlnuoualy  totala  the 
amount  of  electrical  energy  consumed* 

The  printer  componenta  control  the  dlaplay  of  the 
ayatem  outputa*  In  thla  almulatlon  aome  of  the  outputa  are 
printed  at  hourly  lntervala  and  aome  at  dally  lntervala. 

Each  of  the  TRNSYS  componenta  requires  a  aet  of 
defining  parameters  and  Inputs  to  produce  an  output.  The 
parameters  deacrlbe  the  operating  characteristics  of  the 
system  component  being  modeled*  Many  of  these  parameters 
have  been  dlscuased  previously*  A  complete  list  and 
description  of  the  typical  parameters  used  in  this 
simulation  for  each  component  can  be  found  in  Appendix  B. 


CHAPTER  3 


3.1  Presentation  and  Dlaeuaaion  of  Results 

The  collector  efficiency,  which  is  a  Measure  of  the 
useful  solar  energy  being  transferred  to  the  heating  system, 
is  shown  in  Figures  7  and  8  for  various  values  of  collector 
area  and  storage  capecity*  The  collector  efficiency  is 
maximum  where  the  collector  area  is  small  and  the  storage 
capacity  is  large.  This  can  be  shown  also  by  the  collector 
efficiency  equation  when  one  realizes  that  the  inlet  fluid 
temperature  to  the  collector  (from  storage)  will  be  a 
minimum  when  the  ratio  of  collaetor  araa  to  storage  capacity 
is  a  minimum.  In  other  words  looking  at  equation  (1)  again, 

TJ  -  .  69-.  7  [  (Tfi-Ta)/HTJ ,  (1) 

low  values  of  Tfl  (resulting  from  s  small  collector  area  to 
storage  capacity  ratio)  will  reduce  the  second  term  on  the 
right  hand  side  resulting  in  a  larger  value  for  the 
collector  efficiency.  The  collector  efficiency  Increases 
for  Increases  in  the  storage  capacity  and  decreases  in  the 
collector  area  since  both  of  these  trends  decrease  the 
temperature  of  the  fluid  entering  the  collector.  Figure  7 
also  shows  that  increasing  the  storage  capacity  beyond  about 
2  gal/sqft  of  collector  area  is  of  little  benefit  since  the 
collector  efficiency  remains  relatively  constant. 
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FIGURE  8  -  Collector  efficiency  vs  collector  eree  for 
solar  assisted  hast  pump  system  in 
State  College,  Fa* 


The  effect  of  weather  on  the  collector  efficiency  can 
be  seen  by  comparing  the  curves  for  1960  and  1972  in  Figures 


7  and  8*  The  1972,  or  minimum  solar  insolation,  curves 
reflect  higher  collector  efficiencies  for  identical  system 
parameters  (collector  area  and  storage  capacity)*  Thla 
result  la  not  as  easy  to  predict  from  equation  (1)  since  all 
three  variables  on  the  right  side  change.  From  the 
graphical  results,  however,  it  can  be  concluded  that  for 
poor  weather  the  predominant  effect  on  the  efficiency 
equation  is  a  decrease  of  the  second  term  on  the  right  hand 
aide*  This  implies  that  the  magnitude  of  the  drop  in 
entering  fluid  temperature  must  be  larger  than  the  drop  in 
solar  insolation,  indicating  relatively  low  entering  fluid 
temperatures  (storage  tank  temperatures)  would  be  expected 
for  poor  weather  conditions.  Figure  7  also  shows  that 
weather  extremes  have  a  negligible  effect  on  the  collector 
efficiency  of  systems  with  small  collector  area  and  a 
significant  effect  on  the  collector  efficiency  of  systems 
with  large  collector  area. 

Values  of  the  System  Coefficient  of  Performance  (SCOP) 
for  various  collector  area  and  storage  capacity 
configurations  are  shown  in  Figures  9  and  10.  The  SCOP  is 
defined  by  the  following  equation: 

total  heating  load 

SCOP  •  ----- - - — - - - - -  (7) 

electric  energy  used  by  hast  pump  6  aux  heater 
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FIGURE  10  -  System  coefficient  of  performance  vs  collector 
area  for  solar  assisted  heat  pump  system  in 
State  College,  Pa. 
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Following  conventional  practice,  tho  energy  required  for 
fene  end  puaipe  le  eoneidered  negligible  in  coaparlaon  with 
ocher  energy  inputs.  Since  Che  total  heating  load 
requireaenc  for  the  realdence  la  independent  of  the  heating 
systea  configuration,  the  SCOP  will  depend  only  upon  the 
eaount  of  energy  ueed  to  run  the  heat  puap  and  the  energy 
required  for  the  auxiliary  heater*  The  required  energy 
input  froa  theae  two  soureea  ie  reduced,  however,  when  there 
ie  a  contribution  froa  direct  heating. 


In  Figurea  9  and  10  it  can  be  seen  that  all  the  curves 
spproach  s  SCOP  of  3.0  ss  the  oversll  systea  site  (collector 
ares  and  storage  capacity)  inereaaes.  This  is  because  es 
the  storage  capacity  to  collector  area  ratio  lnereasea  the 
aystea  becomes  aore  and  aore  dspendent  upon  the  heat  puap 
alone  (l.e.  the  direct  heating  contribution  and  the 
auxlllery  heat  contribution  to  the  heetlng  load  decrease). 
Specifically,  the  curves  converge  to  s  value  of  3*0  since 
that  la  the  average  heating  COP  of  the  AAF  heat  puap  over 
its  range  of  evaporator  inlet  (storage  tank)  teaperatures. 


Figures  9 

and  10 

both  show 

that 

the 

aaxlaua 

SCOP  is 

obtained  with 

a  large 

collector 

area 

and 

a 

anal  1 

storage 

capacity.  The 

factor 

responslb le 

for 

driving 

the  SCOP  past 

3.0  (the  average  heating  COP  of  the  heat  puap)  is  the  direct 
heating  contribution.  Without  a  direct  heating  contribution 
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the  SCOP  will  be  Halted  to  a  maximum  value  equal  to  the 
heating  COP  of  the  heat  punp.  Also  froa  Figure  10  It  la 
aeen  that  for  a  SCOP  greater  than  3*0  the  aolar  aaeleted 
heat  puap  ayetea  aust  have  aa  a  alnlaua  a  collector  area 
greater  than  about  6S0  sqft. 

The  collector  area  of  630  aqft  In  Figure  10  la  alao 
laportant  In  that  It  la  the  approxlaate  collector  area  for 
which  all  the  different  atorage  capacity  curve*  resulting 
from  1960  weather  data  Intersect  each  other  at  a  SCOP  equal 
to  3.0.  The  significance  of  this  Intersection  Is  that  It 
represents  the  point  at  which  the  total  residence  heat 
requlreaent  is  provided  by  the  heat  puap  and  systea 
perforaance  la  Independent  of  the  storage  capacity.  This 
aeans  that  at  630  sqft  of  collector  area  the  energy  Input 
froa  the  solar  collector  array  la  equal  to  the  energy 
withdrawn  froa  the  storage  tank  by  the  heat  puap,  resulting 
In  neither  the  benefit  of  a  direct  heating  contribution  or 
the  requlreaent  for  auxiliary  heat.  By  extrapolation  it  Is 
estlaated  that  a  aiaillar  intersection  point  will  occur  for 
the  atorage  capacity  curves  resulting  froa  1972  weather  data 
at  about  1800  sqft.  Since  1960  and  1972  represent  the 
extreaes  of  solar  Insolation  It  Is  logical  to  expect  that 
for  the  average  year  the  lnteraectlon  point  will  be 
soaewhere  between  630  and  1800  sqft*  say  1200  sqft.  Thus, 
based  on  a  year  of  average  weather  the  solar  assisted  heat 
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pump  system  for  Staca  Collega  would  raqulrc  a  collector  area 
of  1200  sqft  or  greater  to  obtain  a  SCOP  greater  than  3.0 
(the  heat  pump  heating  COP). 


For  collector  arena  below  the  interaectlon  point  the 
realdenee  heat  load  la  aupplied  essentially  by  the  heat  pump 
and  auxiliary  heat  and  the  SCO P  increases  as  the  storage 
capacity  lncreaaes.  For  collector  arena  above  the 
Interaectlon  point  the  residence  heat  load  la  supplied  by 
the  heat  pump*  auxiliary  heat  and  direct  heat  and  the  SCOP 
decreases  as  the  storage  capacity  Increases. 

This  situation  suggests  thst  a  two  tank  storage  system 
would  be  beneficial.  A  smaller  high  temperature  tank  would 
be  charged  first  by  the  solar  collectors  and  be  used  for 
direct  heating.  A  second  larger  tank  In  the  medium 
temperature  range  would  be  charged  with  the  excess  solar 
energy  and  be  uaed  as  a  heat  source  for  the  heat  pump.  This 
combination  would  allow  the  direct  heating  contribution  to 
be  maximized.  Whether  or  not  the  advantages  of  a  two  tank 
system  outweigh  the  added  coat  and  complexity  would  have  to 
be  studied  in  more  detail.  Figures  9  and  10  also  Indicate 
that  for  a  single  storage  tank  system  the  SCOP  Is  more 
strongly  influenced  by  the  collector  area  than  by  the 
storage  capacity,  since  in  general  the  SCOP  Increases  more 
for  an  Increase  In  collector  area  chan  for  an  Increase  in 
storage  capacity. 
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Comparing  Che  curves  resulting  from  1960  weather  data 
(good)  and  1972  weather  data  (poor)  In  Figures  9  and  10 
shows  that  the  weather  data  used  for  a  solar  assisted  heat 
pump  system  simulation  has  a  major  effect  on  the  resulting 
system  performance  characteristics*  The  SCOP  obtained  for  a 
given  collector  area-storage  capacity  configuration  can  vary 
by  more  than  2502*  If  non-representative  weather  data  are 
used  in  the  design  calculations ,  a  system  design  could 
result  which  is  Inadequate  for  one  situation  or  overdesigned 
for  another. 

Figures  11,  12  and  13  show  tha  precentage  contribution 
made  by  the  heat  pump,  direct  heat  and  auxiliary  heat  to  the 
total  heating  load  for  a  range  of  collector  areas  and 
storage  capacities.  For  a  given  collector  area  the  heat 
pump  contribution  increases  as  the  storage  capacity 
increases.  Accordingly,  the  percentage  of  the  total  heat 
load  supplied  by  auxiliary  heat  decreases.  Looking  at  the 
bar  graphs  for  1960  and  1972  (representing  weather  extremes) 
it  is  seen  that  poor  weather,  as  intuitively  expected. 
Increases  the  auxiliary  heat  requirements.  The  amount  of 
Increase  is  most  severe  around  collector  areas  of  600  sqft 
where  it  averages  about  40Z.  The  percentage  contribution 
made  by  the  heat  pump,  for  a  fixed  storage  capacity  to 
collector  area  ratio,  peaks  in  the  middle  range  of  collector 
areas.  As  predicted  previously  by  Figures  9  and  10,  the 
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FIGURE  11  -  Source  of  heat  for  solar  assisted  heat 
pump  system  In  State  College,  Pa. 
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FIGURE  12  -  Source  of  heat  for  solar  assisted  heat 
pump  system  In  State  College,  Pa, 
at  600  sqft  of  collector  area 
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heat  puap  contribution  la  near  1  00  X  for  a  collector  area 
near  650  sqft  (for  1960  weather  data)*  At  large  collector 
areas  (1000  sqft.  Figure  13)  the  heat  pump  contribution  Is 
smaller  for  good  weather  than  for  poor  weather  because  the 
direct  heat  contribution  lessens  the  heat  requirement  to  be 
provided  by  the  heat  pump* 

Direct  heating  becomes  a  significant  contribution  when 
the  collector  area  Is  large  and  the  storage  capacity  la 
small  because  the  large  amount  of  collected  solar  energy 
only  has  to  heat  a  small  amount  of  water,  resulting  In  a 
high  temperature  level  In  the  storage  tank.  Decreasing  the 
collector  area  causes  the  percentage  contribution  of  direct 
heating  to  decrease  rapidly  so  that  for  small  collector 
areas  the  direct  heating  contribution  is  non-existent*  For 
a  sufficiently  large  collector  area,  increasing  the  storage 
capacity  will  decrease  the  direct  heating  contribution. 
Also,  poor  weather  will  decrease  the  likelihood  of  sny 
direct  hasting  contribution.  In  Figure  13,  for  a  storage 
capacity  of  1  gal/sqft  of  collector,  the  direct  heat 
contribution  for  poor  weather  (1972)  Is  only  SX  compared  to 
46Z  for  good  weather  (1960).  As  Indicated  by  previous 
figures,  the  bar  graphs  of  Figures  11,  12  and  13  show  that 
the  best  performimg  configuration  of  the  solar  assisted  heat 
pump  system  Is  a  large  collector  area  with  a  small  storage 
capacity,  mainly  because  of  ths  advantage  of  more  direct 
best  lng . 
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Table  2  lists  the  energy  usage  of  a  system  with  1000 
aqft  of  collector  area  and  a  storage  capacity  of  1  gal/sqfc 
of  collector.  Energy  data  la  given  for  both  "good"  and 
"poor"  weather.  It  la  clear  from  the  percent  deviations 
Hated  that  weather  has  a  major  Impact  on  the  system 
performance  and  energy  requirements.  Note  that  when  the 
direct  heating  contribution  falls  off  during  poor  weather 
the  majority  of  the  difference  la  made  up  by  auxiliary  heat. 
Several  comparisons  can  be  made  from  Table  2  besides  the 
deviations  shown.  For  Instance,  assuming  an  electric  rate 
of  $.035/kwhr  for  State  College  and  the  system  configuration 
given  in  Table  2,  it  would  coat  between  $28  and  $73 
(depending  upon  the  weather)  to  heat  the  house  for  the  month 
of  December. 


TABLE  2  -  A  one  month  comparison  of  system  performance 
characteristics  for  a  solar  assisted  heat  pump 
system  In  State  College,  Pa.  ;  Collector  area  - 
1000  sqft  and  water  storage  capacity  -  1000  gal 


Dec  1960 

Dec  1972 

Percent 

Weather 

Weather 

Deviation 

Data 

Data 

From  1960 

Collector  Efficiency  (X) 

36 

44 

22.2 

Sys  Coefficient  of  Performance 

5. 14 

1.90 

-6  3.0 

Total  Energy  Collected  (kwhr) 

3648 

1778 

-51.2 

Total  Heating  Load  (kwhr) 

4161 

3950 

-5.1 

Heat  Pump  Contribution  (kwhr) 

2  131 

2308 

16.6 

Direct  Heat  Contribution (kwhr) 

1902 

193 

-89.9 

Aux  Heat  Contribution  (kwhr) 

100 

1242 

1142 

Energy  Input  to  Heat  Pump  (kwhr) 

70  9 

84  1 

18.6 

Energy  Input  to  Aux  Heat  (kwhr) 

100 

1242 

1142 
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Sine*  the  collector  array  la  aounted  on  the  roof  of  the 
residence,  In  reality  an  upper  Halt  of  collector  array 
else  la  aet  by  the  area  of  the  roof.  Taking  thla  Into 
eona Ida  rat  Ion  a  collector  area  of  600  sqft  vaa  ueed  for  the 
paraaetrlc  atudlea  of  room  theraoatat  netting  and  ethylene 
glycol/water  heat  exchanger  alze.  Figure  14  ahova  for  a 
collector  area  of  600  aqft  the  effect  on  SCOP  of  lowering 
the  room  teaperature  to  65  degreea  F  as  the  storage  capacity 
Is  varied.  It  la  Interesting  to  note  that  the  SCOP  la 
Increased  when  the  rooa  teaperature  Is  aalntalned  at  65 
degrees  F  versus  70  degrees  F.  Recalling  that  the  SCOP  Is 
defined  by  equation  (7)  as 

total  heating  load 

SCOP  -  - - - -  (7) 

electric  energy  used  by  heat  pump  6  tux  heater 

one  might  expect  that  since  the  heating  load  la  reduced  by 
about  13Z  the  SCOP  would  decrease,  but  the  reduction  In 
energy  Input  to  the  systea  Is  even  greater.  The  heat  pump 
Input  energy  Is  reduced  by  1 3Z  and  the  auxiliary  heat  Input 
is  reduced  to  nearly  zero  resulting  In  a  net  Increase  In  the 
SCOP.  As  shown  by  the  nearly  horizontal  curves  of  Figure 
14,  and  as  previously  stated  for  a  collector  area  of  600 
aqft,  the  storage  capacity  size  has  little  effect  on  the 
SCOP. 

Table  3  shows  how  lowering  the  rooa  teaperature  to  65 
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degrees  F  effects  Che  energy  inputs  of  the  system*  A 
storsge  cspscity  of  6000  gsllons  vis  chosen  since  Figure  10 
shows  thst  for  600  aqft  of  collector  eres  the  sexiaua  SCOP 
is  obtained  with  s  storsge  cspscity  of  10  gsl/aqft  of 
collector  area.  The  economics  of  lowering  the  room 
thermostat  setting  can  be  obtained  from  the  values  in  Table 
3.  Again  assuming  an  electric  rate  of  $.035/kwhr  for  State 
College  it  is  seen  that  for  the  system  configuration  given 
in  Table  3  approximately  $8  to  $18  can  be  saved  (depending 
upon  the  weather)  for  the  month  of  December  by  lowering  the 
thermostat  setting  S  degrees* 

Tables  4  and  S  are  the  system  performance 
characteristics  for  good  (1960)  and  poor  (1972)  weather 
resulting  from  variations  In  the  ethylene  giycol/wacer  heat 
exchanger  effectiveness.  If  the  capacity  rate  ratio  is  held 
constant  at  a  value  of  1.0  (l.e.  Cmln/Cmax  -  1.0)  the  heat 
exchanger  effectiveness  (£)  is  defined  by  the  following 
equation : 

E  -  NTU/ ( 1+NTU )  (8) 

where  NTU  »  Number  of  transfer  units 

(slse  parameter) 

Thus  by  holding  the  capacity  rate  ratio  constant  at  a  value 
of  1.0,  variations  In  the  heat  sxchanger  effectiveness  (E) 
actually  represent  variations  in  the  slse  of  the  heat 
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TABLE  4  -  System  performance  factors  of  solar  asslated  heat 
pump  system  for  vsrying  values  of  heat  exchanger 
ef fectlvenesa  baaed  on  DEC  1960  weather  data  for 
State  College,  Pa. 


STORAGE  COLLECTOR  CONTRIBUTION  TO  HEAT  LOAD  (X) 


HXCR  TANK 

EFF  TEMP  (  F) 

EFFICIENCY  SYSTEM 
(X)  COP 

HEAT 

PUMP 

DIRECT 

HEAT 

AUX 

HEAT 

(Collector 

Area-600sqft  Storage  Capaclty-600gal) 

.6  115.8 

44.3 

2.68 

73 

14 

13 

.7  117.3 

44.7 

2.71 

72 

15 

13 

.8  118.3 

45.0 

2.74 

71 

16 

13 

(Collector 

Area-600sqft  Storage  Capacity-2 400ga 1) 

.6  87.1 

47.7 

2.93 

99 

0 

1 

.7  88.8 

48.0 

2.94 

99 

0 

1 

.8  90.0 

48.3 

2 . 9  S 

99 

0 

1 

(Collector 

Area-600aqft  Storage  Capaclty-6000gal) 

.6  76.7 

48.5 

2.95 

99 

0 

1 

.7  78.0 

49.0 

2.96 

99 

0 

1 

.8  78.9 

49.4 

2.96 

99 

0 

1 

TABLE  5  -  System  performance  factors  of  solar  asslated  heat 
pump  system  for  varying  values  of  heat  exchanger 
effectiveness  based  on  DEC  1972  weather  data  for 
State  College,  Pa. 


STORAGE  COLLECTOR  CONTRIBUTION  TO  HEAT  LOAD  (X) 

HXGR  TANK  EFFICIENCY  SYSTEM  HEAT  DIRECT  AUX 

EFF  TEMP (F)  (X)  COP  PUMP  HEAT  HEAT 


(Collector  Area-600sqft  Storage  Capaciey-600gal ) 


.6 

40.2 

48.8 

1.47 

46 

1 

53 

.7 

39.8 

49.4 

1.48 

46 

1 

53 

.8 

40.1 

49.8 

1.48 

47 

1 

52 

(Collector 

Area-600sqft  Storage 

Capacity 

-2400gal) 

.6 

40.0 

52.6 

1.62 

57 

0 

43 

.7 

40.0 

53.3 

1.63 

58 

0 

42 

.8 

40.0 

53.8 

1.63 

58 

0 

42 

(Collector 

Area-600sqft  Storage 

Capacity 

■6000ga 1) 

.6 

40.0 

52.0 

1.81 

67 

0 

33 

.7 

40.0 

52.7 

1.82 

68 

0 

32 

.8 

40.0 

53.3 

1.83 

68 

0 

32 
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exchanger.  The  saeller  Che  value  of  E  che  saaller  la  che 
alee  of  che  heac  exchanger.  Ae  can  he  eeen  from  Teblee  4 
end  S  varying  che  heac  exchanger  ef f ace Iveneea  (l.e.  alee) 
of  Che  echylene  glycol/vacer  heeC  exchanger  hee  negligible 
efface  on  che  performance  of  che  syscen.  The  SCOP  lnereasee 
very  ellghcly  ea  che  heac  exchanger  alee  lncreaaea 
regardless  of  che  veacher  for  che  ef f ecclvenesa  range 
acudled. 
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CHAPTER  4 


4.1  Conclusion* 

The  maximum  collector  efficiencies  are  obtained  for 
system  configurations  with  saall  collector  areas  and  large 
storage  capacities  because  the  storage  tenk  temperature  and 
thus  the  teaperature  of  the  fluid  entering  the  collector  la 
kept  low.  Similarly,  the  collector  efficiency  of  a  given 
systea  is  laproved  during  periods  of  poor  solar  insolation 
due  to  the  reduced  entering  fluid  teaperature.  The 
collector  efficiency  Increases  as  the  storage  capacity 
Increases  and  as  the  collector  area  decreases  since  the 
eoablnatlon  of  these  result  in  a  lower  temperature  of  the 
fluid  entering  the  collector.  Heather  extremes  have  little 
affect  on  the  collector  efficiency  of  systems  with  small 
collector  area,  however,  the  effect  on  the  collector 
efficiency  is  significant  for  systems  with  large  collector 
erea.  The  study  results  also  show  that  for  the  solar 
assisted  heat  pump  system  a  storage  capacity  beyond  2 
gal/sqft  of  collector  area  is  of  little  benefit  since  the 
collector  efficiency  remains  relatively  constant. 

For  a  solar  assistad  heat  pump  system  the  SCOP  is 
limited  to  the  heating  COP  of  ths  heat  pump  unless  the 
combination  of  good  solar  insolation  and  adequate  collector 


art*  rtiult  la  a  heating  contribution  froa  direct  heating. 
The  SCOP  will  be  largeat  for  systems  with  large  collector 
areaa  and  email  etorage  capacities.  There  le  a  unique 
collector  area  at  which  the  SCOP  of  the  aolar  aeeleted  heat 
puap  ayetea  will  equal  the  heating  COP  of  the  heat  puap 
(l.e.  100Z  of  the  total  residence  heating  requirement  le 
provided  by  the  heat  puap)  for  a  wide  range  of  storage 
capacities.  Thla  unique  collector  area  Is  dependent  upon 
the  solar  Insolation  received  and  will  vary  In  value 
according  to  the  weather.  For  the  aolar  assisted  heat  puap 
system  aodeled  for  State  College  In  this  study  It  Is 
estlaated  that  about  1200  aqft  of  collector  area  would  be 
required  to  provide  100Z  heat  puap  heating  under  average 
Oeceaber  weather  conditions. 

The  SCOP  of  the  aolar  assisted  heat  puap  system 
Increases  with  Increasing  collector  area  and  In  general  le 
aore  strongly  influenced  by  the  collector  area  than  by  the 
etorage  capacity.  For  a  systea  with  a  collector  area 
greater  than  that  resulting  In  100X  heat  puap  heating,  the 
residence  heating  load  is  supplied  by  contributions  froa  the 
heat  puap,  auxiliary  heat  and  direct  heat  and  the  SCOP 
decreases  as  the  storage  capacity  Increases.  For  a  systea 
with  a  collector  area  smaller  than  that  resulting  In  100X 
heat  puap  heating,  the  residence  heating  load  is  supplied 
essentially  by  the  heat  puap  and  auxilary  heater  and  the 
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SCOP  lnertnti 


as  Che  icorige  capacity  Increases.  This 
situation  suggests  the  use  of  a  two  tank  storage  system  with 
a  high  temperature  tank  being  used  to  maximize  Che  direct 
heat  contribution  and  a  medium  temperature  tank  for  use  as 
the  heat  pump  heat  source. 

For  a  given  collector  area  the  percentage  contribution 
of  the  heat  pump  increases  as  the  storage  capacity  increases 
whereas  the  percentage  contribution  of  the  auxiliary  heater 
decreases.  Direct  heating  becomes  a  significant 
contribution  to  the  total  heating  load  for  solar  assisted 
heat  pump  systems  with  large  collector  areas  and  small 
storage  capacities  since  the  storage  la  heated  to  a  high 
temperature  level.  Decreasing  the  collector  area  causes  the 
percentage  contribution  of  direct  heating  to  decrease 
rapidly  so  that  for  small  collector  areas  the  direct  heating 
contribution  la  non-existent. 

When  designing  or  simulating  solar  assisted  heat  pump 
systems  It  Is  essential  that  accurate  and  realistic  weather 
data  be  used.  The  weather  data  used  has  a  significant 
effect  on  the  performance  characteristics  of  the  system. 
For  the  system  modeled  in  this  study  the  SCOP  varied  by  as 
much  as  2502  for  some  col lector /storage  configurations  when 
"poor"  versus  "good"  weather  data  were  used.  The  percentage 
contributions  of  the  heat  pump,  auxiliary  heater  and  direct 
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heating  to  the  total  heating  load  are  similarly  affected  as 
was  shown  In  Figure  13.  A  simple  economic  analysis  of  a 
system  in  State  College  with  1000  aqft  of  collector  and  1000 
gallons  of  storage  shows  that  monthly  heating  costs  can  vary 
from  $28  to  $73  depending  upon  the  weather. 

Federal  energy  conservation  guidelines  call  for 
reducing  thermostat  settings  to  65  degrees  F  for  winter 
heating.  This  study  found  that  reducing  the  room 
temperature  from  70  degrees  F  to  65  degrees  F  increased  the 
SCOP  approxlmat ely  52  and  reduced  significantly  the  amount 
of  auxiliary  heat  required.  For  a  system  configuration  of 
600  sqft  of  collector  area  and  6000  gallons  of  storage 
capacity  a  simple  cost  analysis  shows  that  In  State  College 
approximately  $8  to  $18  per  month  (depending  upon  the 
weather)  can  be  saved  by  lowering  the  thermostat  setting  5 
degrees . 

The  performance  of  the  solar  assisted  heat  pump  system 
is  not  significantly  affected  by  the  effectiveness  (l.e. 
else)  of  the  ethylene  glycol/  water  heat  exchanger.  Since 
this  is  the  case,  it  appears  that  the  heat  exchanger  could 
be  of  relatively  simple  design  and  low  cost. 

Vhen  reviewing  the  results  presented  In  this  study  it 
must  be  remembered  that  several  assumptions  have  been  made 


In  modeling  and  simulating  the  operation  of  the  residential 


aolar  aaslsted  heat  pump  system.  The  majority  of  these 
assumptions  are  Inherent  in  the  TRNSYS  components  but  others 
are  externally  controlled.  For  Instance  the  effect  of  wind 
velocity  has  been  assumed  negligible  which  in  a  "real  world" 
installation  could  reduce  the  solar  energy  collected  by  the 
collector.  Also  a  ground  reflectance  value  for  bare  ground 
has  been  assumed.  In  a  "real  world"  installation  the 
reflectance  value  might  be  higher  if  the  ground  was  snow 
covered  resulting  in  an  increase  of  reflected  solar 
radiation  received  by  the  collector.  Overall,  however,  it 
is  expected  that  the  model  performance  will  be  aomewhat  more 
optimistic  than  that  of  a  "real  world"  ins tallatlcn . 


Although  domestic  hot  water  heating  was  not  included  in 
the  solar  assisted  heat  pump  model  of  the  residential 
heating  system,  the  collector  array  sized  for  the  solar 
assisted  heat  pump  system  would  most  likely  provide  all  of 
the  domestic  hot  water  requirements  in  the  summer  months. 
During  the  heating  season  the  system  would  be  able  to 
provide  preheating  of  the  domestic  hot  water  much  of  the 
time* 


At  first  glance  the  TRNSYS  modeling  program  seems  very 
complicated  but  actually  the  building  block  modeling 
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technique  Is  fairly  straightforward 


Individual  ayatem 


components  (l.e.  pumps,  solar  collectors,  storage  tank, 

; 

etc,)  are  described  by  component  models  which  are  connected 
together  to  form  a  complete  system.  It  la  better  to  keep 
the  model  aa  simple  as  possible  since  problems,  such  as 
failure  of  the  differential  equation  solution  techniques  to 
converge  within  the  component  models,  are  likely  to  result 
if  the  model  Is  too  elaborate.  Some  modes  of  the  TRNSYS 
components  use  empirical  equations  which  have  been  developed 
in  SI  units  requiring  the  entire  simulation  to  be  run  in  SI 
units  if  these  component  modes  are  used.  If  these  component 
modes  can  be  avoided  in  the  model  the  more  familiar  English 
units  can  be  used.  The  TRNSYS  simulations  require  hourly 
weather  data  as  input.  The  assembly  and  handling  of  auch 
data  can  become  quite  cumbersome  if  punched  cards  are  being 
used,  especially  if  the  simulation  is  run  for  a  full  year. 
Thia  problem  can  be  alleviated  by  using  a  tape  input.  The 
computer  time  required  for  running  a  full  year  simulation  of 
a  solar  assisted  heat  pump  system  (based  on  this  one  month 
run)  will  be  about  200  seconds  and  cost  approximately  $60  if 
output  is  printed  for  each  hour. 

A. 2  Recommendations 

To  obtain  a  complete  picture  of  the  solar  assisted  heat 
pump  system,  a  detailed  economic  analysis  of  the  system 
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alternative*  should  be  performed.  This  would  determine  the 
most  cost  effective  combination  of  such  things  as  collector 
area  and  storage  capacity.  An  economic  analysis  would  also 
be  an  affective  method  to  compare  the  solar  assisted  heat 
pump  system  with  conventional  heating  systems. 


An  entire  year  or  at  least  a  complete  heating  season  of 
weather  data  should  be  used  when  simulating  the  solar 
assisted  heat  pump  system  rather  than  a  single  month.  This 
will  give  a  more  realistic  picture  of  the  system's  seasonal 
performance  characteristics.  Limiting  the  weather  data  to 
the  heating  season  will  significantly  decrease  the  amount  of 
data  and  reduce  computer  costs.  It  Is  Important  however 
that  the  weather  data  be  carefully  selected  to  provide  a 
realistic  representation  of  the  location  being  studied. 


Future  studies  on  the  solar  assisted  heat  pump  system 
should  explore  the  feasibility  and  operating  characteristics 
of  a  two  tank  storage  system  for  optimizing  the  output  of 
the  solar  assisted  heat  pump  system.  Also  the  results 
presented  In  this  report  and  any  future  simulations  using 
the  TRNSYS  program  should  be  compared  with  other  modeling 
techniques  such  as  those  being  developed  by  the 
Architectural  Engineering  Department  of  the  Pennsylvania 
State  University. 
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BUILDING  CRITERIA 


Wall; 


GUm: 


Poore; 


Celling! 


Indoor  Design  Temperature  ■  70F 
Outdoor  Deeign  Temperature  ■  7 *F 


CONSTRUCTION  DATA 


4"  Face  Brick 

3/4"  Air  Gap  (Mean  Temp.  »  SO.  Temp.  Diff.  ■  10,  S  ■  .82) 
Building  Membrane  (Permeable  Felt) 

1/2"  Plywood  (Douglas  Fir) 

3  1/2"  Fiberglass  Insulation 
1/2"  Gypsum  Wall  Board 

2x4  Stud  Wall  Construction.  16"  O.C. 


Area  is  given  on  plan;  e.g.  A  ■  15 
All  Double  Pane 

(1/4"  Air  Gap,  Wood  Sash,  80%  Glass) 


3'  x  7' 

1  1/2"  Solid  Wood 

Metal  Storm  Door  (Front  Door  only) 


6"  Fiberglass  Insulation 
1/2"  Gypsum  Board 

2x6  Ceiling  Joists,  16"  O.C. 
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U-Values: 

Walls: 


1977  Handbook  of 
Fundamentals  Reference 


Table  No. 


Outside  Air  Film  (IS  mph  wind) 

4"  Face  Brick 
3/4"  Air  Gap 

(Mean  Temp.  *  50,  Temp.Dlff.»  10, 
E  =  .  82) 

Building  Membrane  (Permeable  Felt) 
1/2"  Plywood  (Douglas  Fir) 

3  1/2"  Fiberglass  Insulation 
1/2"  Gypsum  Wall  Board 
Inside  Air  Film 

TOTAL  RESISTANCE 


FINAL  U-VALUE 


LAAL. 


Glass: 


Double  Pane 
(1/4"  Air  Gap) 

Adjustment  for  Wood  Sash,  80%  Glass 
(Ux  adjustment) 

FINAL  U-VALUE 


1  1/2"  Solid  Wood 
with  metal  storm  door 

FINAL  U-VALUE 


22.ZT 


Cellin 


Air  Film  (Still  Air,  E  ■  .  90) 

6"  Fiberglass  Insulation 
1/2"  Gypsum  Wall  Board 
Inside  Air  Film 

TOTAL  RESISTANCE 

FINAL  U-VALUE 

Kitchen  Wall  Against  Garage: 

Air  Film  (Still  Air,  E  •  .  90) 

3  1/2"  Fiberglass  Insulation 
1/2"  Gypsum  Wall  Board 
Inside  Air  Film 

TOTAL  RESISTANCE 

FINAL  U-VALUE 
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.6/  22.// 

11,00  Z2./4 


Z3.4>7 


,  //itf.g  - 

-  -  -  fcv- 
/7.g/ 

.at 


COOLING  LOAD  SUMMARY 


Walls 

Windows 

Doors 

Ceiling 

Floor 

Occupants 

Cooking  . 

WAIL  Infiltration 

Total 


/g£ 

52 Lk~ 
Jl?=L- 

'6l7L, 
_ o 

700 

’ zoo 

VJ± 

&c4 


Btuh 


Btuh,  Total  Sensible  Heat  (TS 


Total  Heat  •  1. 3  (TSH)  • _ 


Equipment  Size 


_  Total  Heat.  Btuh  _  /  /L/i- 

•  n  "A AX  B -T”  TTT--  —  ■  _ /'IT. 


tons 
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APPENDIX  B 


DEFINITION  OF  TERMS  USED  IN  TRNSYS  COMPUTER  MODEL  FLOW  DIAGRAM 

DELTAU  -  Total  change  in  storage  tank  Internal  energy 

H  -  Solar  radiation  (total)  on  a  horizontal  aurface 

Hb  -  Solar  radiation  (beam)  on  a  horizontal  aurface 

Hd  -  Solar  radiation  (diffuae)  on  a  horizontal  aurface 

HT  -  Solar  radiation  (total)  on  a  tilted  aurface 

RB  -  Solar  radiation  (beam)  on  a  tilted  aurface 

HD  -  Solar  radiation  (diffuae)  on  a  tilted  aurface 

HI  -  Fraction  of  tlmestep  that  heat  pump  operated  in 

water-air  heating  mode 

H2  -  Fraction  of  tlmestep  that  haat  pump  operated  in 

alr-alr  heating  mode 

H3  -  Fraction  of  tlmestep  that  heat  pump  operated  in 

cooling  mode 

IGAM  -  Heat  pump  master  control  Integer;  O-off,  1-on 

Ml  -  Flowrate  of  water  entering  heat  pump  evaporator 

Mo  -  Flowrate  of  water  returning  from  heat  pump 

evaporator 

MODE  -  Mode  of  heat  pump  operation 

OPHRS  -  Total  time  heat  pump  operated  in  heating  mode 

QA1  -  Instantaneous  energy  absorbed  from  liquid  source  in 

water-air  heating  mode 

QABS  -  Total  energy  absorbed  by  heat  pump 

QAC  -  Instantaneous  energy  removed  from  load  in  cooling 
mode 

QAUXH  -  Auxiliary  energy  delivered  to  room  in  heating  mode 

QDH  -  Energy  delivered  to  room  in  direct  heating  mode 

QDUMP  -  Rate  at  which  energy  la  discarded  by  relief  valve 

QENV  -  Rate  of  energy  lose  from  storage  tank  to 

surroundings 


QL  -  In s c an t ane oua  rooa  heating  load 

QLOAD  -  Hourly  heating  load 

QR1  -  Instantaneous  energy  delivered  to  load  in  vatcr-elr 
heating  mode 

QR2  -  Instantaneous  energy  delivered  to  load  in  air-air 
heating  node 

QRC  -  Energy  rejected  to  sink  in  cooling  mode 
QREJ  -  Total  energy  delivered  by  heat  pump 

QSHRTC  -  Cooling  capacity  required  but  not  met  by  heat  pump 

QTANK  -  Rate  of  energy  delivery  from  storage  tank  to  heat 
pump 

Qu  -  Rate  at  which  energy  is  transferred  to  storage  tank 

QUSE  -  Hourly  useful  energy  collected 

QU  -  Total  work  input  to  heat  pump 

RAD  -  Hourly  incident  solar  radiation  on  collector 

T  -  Temperature  of  storage  tank 

Ta  -  Ambient  temperature 

TENV  -  Temperature  of  ground  surrounding  storage  tank 

T1  -  Temperature  of  water  entering  heat  pump  evaporator 

To  -  Temperature  of  water  leaving  heat  pump  evaporator 

TOTQL  -  Total  heating  load 

TOTQU  -  Total  useful  energy  collected 

TOTRAD  -  Total  incident  solar  radiation  on  collector 

TSTOR  -  Hourly  storage  tank  temperature 

WAC  -  Instantaneous  work  input  to  heat  pump  in  cooling 
mode 

UAH  -  Instantaneous  work  input  to  heat  pump  in  heating 
mode 

WORKIN  -  Total  work  input  to  heat  pump  in  heating  mode 

AE  -  Change  in  Internal  energy  of  water  In  storage  tank 
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-  Angle  of  incidence  of  eoler  redietion  on  a  tilted 
surface 


1 

1 


j 
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TYPICAL  TRNSYS  COMPONENT  PARAMETERS 


(A)  UNIT  1,  TYPE  9  READER 


Parameters 


Descript  ion 


Value 


1 

2 

3 

A 

5 


N 

At 

1 

Mi 
A 1 


number  of  values  to  be  read 
time  interval  at  which  data  is 
provided 

position  of  value  on  card  for 
which  units  are  to  be  converted 
multiplication  factor  for  that  value 
addition  factor  for  that  value 


2 

1  hr 

2 

3.686 

0 


(parameters  3.4&S  convert  Langleys  to  Btu/hr-sqft) 


6  L  -  logical  unit  number  of  data  file  8 


(B)  UNIT  2,  TYPE  16  RADIATION  PROCESSOR,  MODE  3 


Parameter 


Descr ipt ion 


Value 


1 

2 

3 

4 

5 

6 
7 


n 

4 

a 

i 

Sc 

P 


mode 

day  of  year  at  start  of  simulation 
lat ltude 

slope  of  collector  surface 
orientation  angle  of  collector 
solar  constant 
ground  reflectance 


3 

335 

40.5  deg 

50.5  deg 
0  deg 

429  Btu/hr-sqft 

.2 


(C)  UNIT  3,  TYPE  21  LIQUID  C0LLECT0R-ST0RACE  MODULE 


Parameter 


Descr lption 


Value 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

1 1 
12 

13 

14 


A  -  collector  area 

Me  -  collector  fluid  flowrate 

Cpc-  collector  fluid  specific  heat 

Ah  -  heat  exchanger  to  tank  flowrate 

Cph-  specific  heat  of  fluid  in 
storage  tank 

F'  -  collector  efficiency  factor 
Tb<  -  transmit  tance-absorptance 
pr oduc  t 

UL  -  collector  loss  coefficient 
Tmax-  temperature  at  which  tank 
relief  valve  opens 
6  -  collect or-tank  heat  exchanger 

effect ivenes  s 

V  -  volume  of  storage  tank 
P  -  tank  fluid  density 
U  -  tank  energy  loss  coefficient 
r  -  ratio  of  tank  height  to  tank 
diameter 


600  a  qf  t 
7500  lb/hr 
.792  Btu/lb-F 
5940  lb/hr 
1.0  Btu/lb-F 


.99 

.684 

.7  Btu /hr-aqf t-F 

2  12  deg  F 

.7 

80.1  cuf  t 

62 . 4  lb/cuf  t 

.05  Btu /hr-sqf t-F 

3 
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IS  To  -  initial  tank  fluid  tamperatura  70  deg  F 


(D)  UNIT  4,  TYPE  20  WATER-TO-AIR  HEAT  PUMP 


Parameter 


Description 


Value 


1 

CPI  -  specific  heat  of  liquid 

entering  heat  exchanger  A 

1.0  Btu/lb-F 

2 

MD0T1  -  mass  flow  rate  of  liquid 
entering  heat  exchanger  A 

2500  lb/hr 

3 

TROOMH  -  constant  room  temperature 
for  heating  mode 

70  deg  F 

4 

TROOMC  -  constant  room  temperature 
for  cooling  mode 

77  deg  F 

5 

TMIN1  -  minimum  liquid  source  temp 
for  water-to*air  heat  pump 

40  deg  F 

6 

TM1N2  -  minimum  air  source  temp  for 
air-to-air  heat  pump 

1000  deg  F 

7 

NDATAH  -  number  of  equally  spaced 
heating  data  points 

8 

8 

NDATAC  -  number  of  equally  spaced 
cooling  data  points 

8 

9 

ISAVEH  -  logical  unit  number  of 
heating  data  file 

31 

10 

ISAVEC  -  logical  unit  number  of 
cooling  data  file 

32 

11 

TCOOL  -  minimum  ambient  temp  vhen 
cooling  is  allowed 

1000  deg  F 

12 

CAIR  -  MCp  product  for  room  air 
flow  thru  heat  exchangers 

B,C  and  D 

648  Btu/hr-F 

13 

EFF  -  effectiveness  of  heat 

exchanger  B 

.7 

14 

TSET  -  minimum  liquid  source  temp 
for  direct  heating  from 
liquid  source 

110  deg  F 

15 

ICOOL  -  cooling  condenser  selection 
integer,  0-heat  exchanger  A 

0 

16 

THEAT  -  maximum  ambient  temp  vhen 
heating  is  allowed 

65  deg  F 

UNIT  5, 

TYPE  12  HEATING  LOAD,  MODE  3 

Parameter  Description 

Value 

1 

mode 

3 

2 

UA  -  heating  requirement  of  house 
(energy /degree  hour) 

-511  Btu/hr-F 

3 

TR  -  room  temperature 

70  deg  F 

60 


(F)  UNIT  6,  TYPE  24  QUANTITY  INTEGRATOR 
No  parameter*  required 

(C)  UNIT  7,  TYPE  25  PRINTER  (HOURLY) 

Parameter  Description 

I  Tp  -  time  interval  at  which  printing 

ia  to  occur 


Value 


1  hr 


(H)  UNIT  8,  TYPE  25  PRINTER  (24  hour) 
Parameter  Deacr lotion 

1 


Tp  -  time  interval  at  which  printing 
is  to  occur 


Value 
24  hr 
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T4NSVS  -  *  IHNSIINI  SINUltTION  ptn'.P  t>* 

/•UN  TMt  101.44  6MMC.V  l*S  4T  f HI  JNI1C4SII*  OF  ■  I  }(  UN  VJN _ 

K4SI0N  S.t  1/2  7/76 


Sl'xUltl  ION 
iixe«»'.CEs 


T.*Mt  o*  J.IOOI-M 


-1.0006-02  -1 .0001-01 


UNIT  I  17*6  1  Ft 

MIU1EICIS  s 
2.1906  III  l.CJO£  10 


Pfe4THE4  04)4  *640(4 


2.0306  90 


1.6*6*  00 


9  l»»E  16  44 

nninm  _7 
1.0006  30  '  I.110E  02 

i'i»or$  I 
I.  9 
0.0 


6601  6? I  ON  P40C6SSQ4  i  loot  1 


4.0106  01 


1.010k  01 


.0.3 _ 0.0006  00 


6. 30 OE  02 


2.0)0f-0l 


1  TT*f  21  II 
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2.  1201  02  T.  )91E  •  01 
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2.  I  I.  I 
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